Orbital polarons versus itinerant Cg electrons in doped manganites 
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We study an effective one-dimensional (ID) orbital t-J modei derived for strongly correlated Cg 
electrons in doped manganites. The ferromagnetic spin order at half filling is supported by orbital 
superexchange oc J which stabilizes orbital order with alternating — y^ and Sz^ — orbitals. 
In a doped system it competes with the kinetic energy oc t. When a single hole is doped to a 
half-filled chain, its motion is hindered and a localized orbital polaron is formed. An increasing 
doping generates either separated polarons or phase separation into hole-rich and hole-poor regions, 
and eventually polarizes the orbitals and gives a metallic phase with occupied 3z^ — r^ orbitals. 
This crossover, investigated by exact diagonalization at zero temperature, is demonstrated both 
by the behavior of correlation functions and by spectral properties, showing that the orbital chain 
with Ising superexchange is more classical and thus radically different from the ID spin t-J model. 
At finite temperature we derive and investigate an effective ID orbital model using a combination 
of exact diagonalization with classical Monte-Carlo for spin correlations. A competition between 
the antiferromagnetic and ferromagnetic spin order was established at half filling, and localized 
polarons were found for antiferromagnetic interactions at low hole doping. Finally, we clarify that 
the Jahn- Teller alternating potential stabilizes the orbital order with staggered orbitals, inducing the 
ferromagnetic spin order and enhancing the localized features in the excitation spectra. Implications 
of these findings for colossal magnetoresistance manganites are discussed. 
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I. INTRODUCTION 

Doped perovskite manganese oxides Ri_a;Aj;Mn03, 
where R and A are rare earth and alkaline earth ions, 
have attracted increasing attention because they show a 
rich variety of electronic, magnetic and structural phe- 
nomena, and several different types of ordered phasessA- 
To explain the colossal magnetoresistence^ (CMR) and 
metal-insulator transition observed in these compounds 
as a function of either doping x or temperature T, which 
suggests their potential technological applications, one 
has to go beyond the simple double exchange model 
of ZenerfS- and investigate a complex interplay between 
magnetic, orbital, and lattice degrees of freedom, as well 
as the conditions for the itinerant behavior of strongly 
correlated Cg electrons. 

The theoretical challenge is to understand the prop- 
erties of doped manganites in terms of the dynamics of 
correlated eg electrons which involves their orbital de- 
grees of freedomii Although several features such as 
the transition to ferromagnetic (FM) order under dop- 
ing, and the magnetic excitations in the FM phase, were 
qualitatively reproduced within the Kondo model which 
assumes a nondegenerate conduction bandfii^iSiLSiSiiSiii 
it is evident that the degeneracy of Cg orbitals plays a ma- 
jor role in the double exchange model, in the transport 
properties of doped manganites, and in the CMR effect 
itselfii2iiii4 In the undoped LaMnOa charge fluctuations 
are suppressed by large on-site Coulomb interaction U, 
leading to superexchange which involves both spin and 



orbital eg degrees of freedom. 

Purely electronic superexchange models describing 
the orbital order of eg electrons were suggested early 
on by Kugel and Khomskii,.^^ and are of great inter- 
est recentlyii^ In cuprates, like in KCuFs, they lead 
to enhanced quantum effects^ while in manganites 
these models, which originate from several charge ex- 
citations due to either eg or t2g electrons, are richer 
and include interactions between large and more clas- 
sical spins ii2ii2i22i2i These interactions, together with 
the Jahn- Teller (JT) effect, stabilize A-type antiferro- 
magnetic (AF) phase coexisting with orbital order in 
undoped LaMnOa. While the shape of the occupied 
orbitals in the low-temperature phase of LaMnOs is 
still controversial, the two-sublattice orbital alternation 
in an orbital ordered state within the FM planes of 
LaMnOa is well establishediSS When the spin dynam- 
ics in the FM planes of LaMnOa is frozen at low tem- 
perature, one can constrain an effective model to purely 
orbital superexchang^^ instead of considering more com- 
plex complete spin-orbital model^ The superexchange 
orbital interactions favor then alternating directional 
(3z^ — r^-like) and planar (a;^ — y^-like) orbitals along 
every cubic directioni^ 

A considerable simplification is allowed in the FM 
phase at finite doping and at T = 0, where the spins 
are aligned and an effective t- J-like charge-orbital model 
is sufficient. Experimentally, at doping higher than 
X ~ 0.10 one finds a FM insulating phase, followed by 
a metallic phase. The microscopic reasons of this be- 
havior are intriguing — it contradicts the usual spin po- 
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laronic picture of a FM phase. A puzzling competition 
between the insulating and metallic behavior within the 
FM phase was also reported for Lao.88Sro.i2Mn03iS^ 

Here we will concentrate on the generic behavior due 
to Cg orbital degrees of freedom in a FM phase, and 
we will treat the spin dynamics classically, as usually 
done in the double exchange model. '^■^'^ This approach 
is complementary to focusing on the quantum effects in 
double exchange, which was presented recentlyi2& Our 
aim is to study the correlation functions and spectral 
properties of a one-dimensional (ID) orbital t-J model 
by exact diagonalization at zero temperature {T = 0), 
and by classical Monte Carlo simulations of spin corre- 
lations at finite temperature in order to establish the 
consequences of orbital order and orbital dynamics in 
doped manganites. In the strongly correlated regime the 
charge dynamics couples to orbital excitations and can 
be described in terms of the orbital t-J model for doped 
manganites. As in the spin t-J model, a doped 
hole moves in an orbital ordered state by dressing itself 
with orbital excitationsi^l However, the structure of the 
quasiparticle (orbital polaron^?) is here quite different 
from that derived from the t-J modeljSS as the orbital 
model is more classical. 

As a reference system we use a ID spin t-J model for 
which the behavior of a single hole is well understood. 
A hole created in a ID Neel state is mobile and may be 
thought as decaying into a magnetic domain wall. In 
fact, it gives a charged domain wall already after a sin- 
gle hop and leaves behind a solitonic defect from which 
it separates and next propagates independently. This is 
the simplest visualization of hole-spin separation in a ID 
system. The familar string picture and quasiparticles on 
the energy scale of ~ 2 J are then recovered when a stag- 
gered magnetic field is applied and suppresses domains 
of reversed spins. '^''■^^ 

We show below that an analogous phenomenon of hole- 
orbital separation does not occur for the orbital degrees 
of freedom, but instead a single hole is trapped in a ID 
chain when orbitals alternate. This explains why an insu- 
lating behavior may extend to finite doping. A staggered 
field has a well established physical origin in this case, 
and could follow from the frozen JT modes with alternat- 
ing oxygen distortions along the chain itself, and on the 
bonds perpendicular to its direction, when the chain is 
embedded within a three-dimensional (3D) crystal. This 
interaction may play also an important role for the trans- 
port properties of lightly doped manganites and is diffi- 
cult to separate from on-site Coulomb interaction U.^- 
Therefore, it was even suggested that models including 
Hund's exchange Jh between conductions electrons and 
t2g spins, but neglecting Coulomb interaction U, could 
capture the essential physics of manganites42i We will 
investigate below the JT term which leads to a staggered 
field acting on the orbitals. Note, however, that the ID 
case is special as static JT distortions, preventing fluc- 
tuations of occupied orbitals, suppress completely hole 
motion along the chain, while in a two-dimensional (2D) 



model they lead instead to an enhanced coherent com- 
ponent in the hole motion mainly due to suppression of 
incoherent processesi^S 

Another feature characteristic of the orbital physics 
and different from spin models is that the lattice responds 
to the doping. As pointed out by Kilian and Khaliullin^^i 
the breathing motion of the MnOg octahedra provides a 
strong tendency towards hole localization in the presence 
of singly occupied (almost) degenerate Cg levels, as shown 
in a 2D modeli^ A static hole polarizes then the orbitals 
occupied by Cg electrons in its neighborhood, and this 
polarization is expected to happen in addition to the ef- 
fects promoted by the orbital superexchange interactions. 
We also consider this interaction in the present ID model 
and show that it simply renormalizes the superexchange 
interaction, and thus it could lead to qualitatively new 
features and be of more importance only in higher di- 
mension. 

The paper is organized as follows. In Sec. II we present 
the effective t-J orbital model for Cg electrons moving 
along a ID chain in a FM plane of, e.g., Lai_a;Sra;Mn03. 
In Sec. Ill this model is analyzed first qualitatively 
by comparing the energies of two phases: (i)insulating 
phase with localized eg polarons, and (ii) metallic phase 
with itinerant carriers. Next we analyze the correlation 
functions and the spectral properties of a finite chain of 
TV = 14 sites filled by up to five holes, and demonstrate 
a crossover from the insulating to metallic phase under 
increasing doping. At finite temperature T we derive an 
effective t-J orbital model, with electron hopping and 
orbital interactions depending on actual configuration of 
core spins (Sec. IV). The interrelation between spin and 
orbital order, the correlations around a doped hole, and 
the evolution of spectral properties under increasing dop- 
ing are next investigated using a combination of exact di- 
agonalization with Monte Carlo simulations. Finally, we 
summarize the results and present general conclusions in 
Sec. V. 

II. ORBITAL t-J MODEL 

We consider the ID orbital t-J model, 

no = Ht + Hj + HjT, (1) 

obtained for Cg electrons in FM manganites at T = 0. 
The Hamiltonian Q acts in the restricted Hilbert space 
without double occupancies. Due to the absence of the 
SU(2) symmetry, the kinetic energy in the eg band takes 
a form which depends on the used orbital basis. For 
the present ID model it is most convenient to consider 
a chain along c axis and to use the usual orbital basis 
{x^ — y'^,3z^ — r^}, for which we introduce a compact 
notation, 

\:c)^^\x^~y^), \z)^^^\3z^-r'). (2) 

A chain along a or 6 axis could also be analyzed using this 
basisp^ but in each case one obtains a simpler and more 
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transparent interpretation of the results with a basis con- 
sisting of a directional orbital along this particular cubic 
axis, and an orthogonal to it planar orbital, e.g. — 
and — for a axis. 

The model given by Eq. |^ stands for a chain com- 
posed of Mn^+ and Mn'*+ ions coupled by the effective 
hopping which originates from the transitions over oxy- 
gens. Due to the symmetry, the |a;} electrons cannot 
move, and the hopping term Ht allows only for intersite 
transitions between a site occupied by a \z) electron and 
an empty neighboring site, 

iit = + A+\,zC^z). (3) 

i 

where an operator c|^^ = c|^(l — n^^:) creates b.\z) electron 
when site i is unoccupied by an \x) electron. The hopping 
Ht describes thus spinless fermions with an orbital flavor 
in a restricted Hilbert space, in analogy to the original 
t-J model in spin spacCj-^^ but only one (z) component is 
here itinerant^ while the other {x) one is immobile and 
hinders the motion of \z) electrons along the chain. 

The superexchange in undoped LaMnOa is given by 
a superposition of several terms which originate from 
charge excitations due to either Cg or t2g electron 
hopping?? — they are discussed in Sec. IIV Al While the 
t2g part is AF, the eg terms favor either FM or AF spin 
order on a bond (ij), depending on the pair of occupied 
eg orbitals at sites i and j. For the realistic parameters of 
LaMnOa, one finds the A- type AF order in the ground 
state, with FM planes staggered along the third direc- 
tion. Taking a cubic c direction within a FM plane, the 
superexchange expression simplifies enormously. ^'^ Treat- 
ing large S* = 2 spins at Mn^"*" ions classically, all AF 
terms drop out at T = 0, and the remaining orbital su- 
perexchange favors alternating eg orbitals. '^^ Therefore, 
the superexchange interactions reduce then to the purely 
orbital interactions which favor alternating directional 
(3z^ — r^-like) and planar (x^ — y^-like) orbitals along 
every cubic direction. ■^'^ In the present ID model one finds 

Hj = 2JY, (t^'T^+i - , (4) 

i 

where operators 

i;^ = iaf = i(n„-n,,), (5) 

stand for orbital pseudospins T = 1/2, with two eigen- 
states defined in Eq. and fii = c\^c^^ -\- cf^Cj^ is an 
electron number operator in the restricted Hilbert space. 
The superexchange constant J = /e{^Ai) is then given 
by the high-spin excitation energy e{^Ai)w& 

The last term in Eq. stands for the staggered field 
induced by the cooperative JT effect, 

i/jT = 2EiT eM^T^Ri)Tt. (6) 

i 

considered also in Ref. l37l and supporting the alter- 
nating orbital order. It gives an energy gain _Ejt per 



site in the ground state of an undoped ID chain, and fol- 
lows from the alternating oxygen distortions around man- 
ganese ions in LaMnOsi^ The present simplified form 
[Eq. ©] of a general expression'^^, which depends on the 
type of oxygen distortion, is sufficient in the ID model. 

At finite doping x = 1 — n, where n is an average eg 
electron number per site, the present t-J model gives an 
interesting problem, with competing tendencies towards 
orbital alternation in insulating state on one hand, de- 
scribed by orbital correlations, 

T{n) = (7) 

and uniform \z) polarization in metallic state on the 
other. The nature of the ground state obtained at finite 
doping is best investigated by considering a few char- 
acteristic correlation functions, containing information 
about the orbital state and about the orbital order be- 
tween two nearest neighbor sites, both at distance n from 
a hole, and about the orbital correlations across a hole. 
To optimize them we introduce: 

P{n) = {n^n+^), (8) 
R{n) = (n.i;!,„7;!,„+i), (9) 
Z = {Tt_,n,Tt+,), (10) 

where fii is a hole number operator at site i, 

'iT'i = ^ ~ '^iz'^iz ~ Ax^ix- (11) 

These correlation functions will be discussed in Sees. 
Ill and IV. Note that the kinematical constraint gives 
P(0) = i?(0) = 0. 

We compare the orbital correlations and the spec- 
tral properties obtained for the orbital chain, with these 
found for a ID t-J spin model for 5 = 1/2 spins, 

i(T 

+ 2 J ^ ■ Si+i - inifij+i^ 

i 

+ 2/i,^exp(i^i?05f, (12) 

i 

where — c\^{l — rii^-^) is a creation operator of an 
electron with spin a at site i in the restricted Hilbert 
space. In contrast to the orbital model, here the elec- 
trons with both spin flavors are mobile and exchange be- 
tween the sites i and i -t- 1, so the interaction is a scalar 
product Si ■ Si+i instead of the Ising term TfTl^^. For 
convenience, we use the same units as in Eqs. (0} and 
©; here J — 2t^ /U with U being the excitation energy. 
The staggered field oc hs simulates the long-range spin 
order present in a 2D model^Si^ and plays a similar role 
to the JT field in the present orbital model. For doped 
spin chain we consider analogous correlation functions to: 
orbital-orbital (0, hole-orbital (|H1), orbital order at dis- 
tance n from the hole 0, and around a hole (|10|l . with 
spin operators Sf in place of [Eq. (jSJ] . 
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In the next Section we report the results obtained by 
Lanczos diagonaUzation of an A'^ = 14 orbital chain with 
periodic boundary conditions and different electron fill- 
ing at T = 0, and compare them to more familiar spin 
physics. We have verified that the results concerning the 
crossover to metallic phase are representative and do not 
depend on the chain length in any significant way. 

III. NUMERICAL RESULTS AT T = 

A. Analytic estimation of the crossover from 
orbital polarons to itinerant electrons 

The ground state of the orbital t-J model ^ depends 
on two parameters: J/t and Ej^/t, and on hole doping 
X = 1 — n. We first investigate the ground state (at 
T = 0) changing these parameters. A realistic value of 
t ~ 0.4 — 0.5 eV was estimated for manganites using the 
charge-transfer modeliSS Taking this as an energy unit, 
and the spectroscopic value for the energy of the high- 
spin excitation e{^Ai) — U — 3Jh — 3.8 eV, this leads 
to J/t ~ 1/8. We will consider also Ejt > 0, which 
promotes localized behavior of eg electrons. 

At half- filling one finds alternating \x)/\z) orbitals in 
the ground state at J > 0, with classical intersite corre- 
lations T{n) = (— 1)"/4. If a single hole is then doped 
to orbitally ordered state at an \x) site, it can delocalize 
within a box consisting of three sites, as each neighboring 
\z) electron can interchange with the hole [see Fig. da)]. 
This state is therefore favored over doping at |z) site, and 
one can easily determine the energies of a doped hole in 
an antibonding and bonding state, 

1 - 1 / - \ 1/2 

Eii,,± ^ -J ± -i^J^ + 8t^) , (13) 

with J ^ J + 2EjT standing for an excitation energy of 
the configuration with a hole moved to a left (right) site 
within a three-site cluster shown in Fig. Q^a). 

As long as the holes may be doped into separated three- 
site units (for x < 0.25), the total energy (per site) of an 
insulating phase follows from the weighted contributions 
of the undoped orbital ordered regions, and doped holes 
occupying the bonding states of individual clusters, with 
energy Ei,,^^ 

Ei = -(1 - x)EjT - (1 - 2x)J + xEih,^. (14) 

The energy of a metallic phase, which contains only itin- 
erant \z) electrons for n > 0.5 [see Fig. ^b)], 

Em = - [ Efc.- + - / ^k.+ , (15) 

Jk<Tr/2 JkKkp 

is obtained by integrating over the occupied one-particle 
band states ek,±- At finite £'jt one finds, 

£fe,± = ±^Ej^+4£^ cos^k. (16) 
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FIG. 1: (color online) Competition between localized holes 
within an alternating orbital (AO) phase, with staggered lo- 
calized ja;) (shadded vertical boxes) and itinerant \z) (empty 
horizontal boxes) orbitals (a), and itinerant holes in ferro or- 
bital (FO) \z) phase (b), in a ID chain along c axis. A hole 
(shadded circle) is confined to a cluster of three sites in the 
AO state, while it is delocalized over the entire chain in the 
FO state, as indicated by dashed boxes. Part (c) shows a 
qualitative phase diagram of the ID orbital t-J model 0, 
with a critical concentration Xc separating the AO and FO 
phases, as obtained for _Ejt = and: A = (solid line), 
A = t (dashed line). Filled (empty) circles in (c) show the 
AO (FO) states found by exact diagonalization of an A'^ = 14 
chain. 



The superexchange cx J does not contribute to Em- 

It is now of interest to compare the energy of two ex- 
treme situations: (i) an insulating phase with holes lo- 
calized within three-site clusters, [Eq. H14I) ]. with (ii) a 
metallic (itinerant) phase [Eq. (|15|l ]. One finds that at 
fixed J and iJjx = 0, an insulator-metal transition takes 
place when the hole concentration x increases. A critical 
concentration for this transition increases with increasing 
J/t [Fig. n^c)]. The data points obtained from the exact 
diagonalization of an = 14 chain agree with the ana- 
lytic estimation at J/t =1/12 and 1/8, while for a larger 
value of J/t = 1/4 the region of the insulating AO phase 
is more extended in the chain. Note, however, that the 
present analytic estimate is anyway only qualitative as 
the energy of the AO phase can be evaluated using a su- 
perposition of holes confined to three-site clusters only up 
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to X = 0.25. An extended region of stability of the AO 
phase results here from larger clusters of the itinerant 
phase which are still separated by immobile \x) states; 
such larger itinerant units occur already for x < 0.25 
at low values of J/t, as the gain in the kinetic energy 
for a hole moving within a larger cluster approaches fast 
the metallic limit with increasing cluster size, and may 
be easily compensated when a few exchange bonds are 
created. Indeed, the data show [see Fig. ^c)] that the 
ground state for x — 1/14 remains insulating with two 
electrons at J = 0.02t, and becomes metallic only for 
J < 0.02i. 

A simple estimation of the metal-insulator transition 
is also possible when a polaronic polarization around a 
hole^'^ is included in Eq. 



Ha = -Ay^ni(n,_i, 



(17) 



where njz = c|^i z^i±i z the electron number operator 
in the restricted space without double occupancies, and 
fii is a hole number operator pi|l . At i?jT = a similar 
expression to Eq. H13|l is then obtained, with J = J-|- A, 
and the range of an insulating AO phase shrinks in the 
phase diagram of Fig. ^c). In fact, the polarization 
around doped holes is optimized in a metallic phase, and 
this result is special to the ID model. Therefore, we shall 
not consider this interaction further, as it gives qualita- 
tively the same results as the orbital t-J model at A = 0, 
but with a somewhat reduced value of J. 

At finite JT energy _Bjt the situation changes in a dras- 
tic way. Even for relatively small Ejt — 0.25t the AO 
phase is stabilized in the entire range of doping shown in 
Fig. IHc). Note that even somewhat higher values could 
be more appropriate for realistic 3D manganitesi^^i^S 
This suggests that the JT effect may indeed play a very 
important role in manganites and stabilize an insulating 
phase with orbital order in a broad regime of doping; we 
address this question in Sec. IIII El 



B. Orbital order at half filling 

The superexchange interaction in the orbital model 
0) is Ising-type, and therefore the orbital order in the 
undoped system at finite J is perfect, with alternat- 
ing occupied \x) and \z) orbitals along the chain, and 
{T^T^j^-y) = —0.25. The classical character of this ground 
state is reflected in the one-hole excitation spectra. If a 
single hole is added at half filling (n = 1), it may be 
doped either at an \x) or at a \z) site. A hole doped at 
a \z) site is immobile and the energy of the final state 
is higher by 2 J than that of the initial state, as two ex- 
change bonds are removed. This excitation appears as a 
localized peak at hole binding energy — 2J in Fig. 
In contrast, when a hole replaces an electron at site 
I, the resulting state is not an eigenstate of TYq, the hole 
delocalizes over a three-site cluster including also occu- 
pied \z) orbitals at neighboring sites i — 1 and i + 1 [Fig. 
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FIG. 2: (color online) Spectral functions A{k, ui) at half filling 
for an A'^ = 14 site chain: (a) orbital model, with hole excita- 
tions within occupied |«)/|a:) orbital shown by solid/dashed 
lines, respectively; (b) spin t-J model. Insets show the den- 
sity of states n{uj). Parameters: J — 0.125t, Ejt — 0, hs — 0. 
We assumed a finite broadening of the peaks hy S = O.lt. 



Wi^)]i and the bonding and antibonding hole states l|13() 
contribute. They give two maxima in the excitation spec- 
trum of Fig. [2Ia), with energies Eih^± — 2 J. As the hole 
cannot hop over the sites occupied by \x) electrons, these 
two states are again localized within a cluster, and lead 
to fc-independent maxima in the spectral function. 

Strictly speaking, three-site terms, similar to those 
known from the derivation of the spin t-J model;^ oc- 
cur as well in the orbital model. One might expect that 
such terms, oc Jcl^^^ ^hi^Cizpi.z, would change the spec- 
tral functions as a hole created at \z) site could then hop 
to its second neighbors and interchange with immobile 
electrons in |a:) orbitals. However, we have verified that 
these processes are of importance only for J ~ but do 
not lead to any significant changes of the spectral prop- 
erties in the physically interesting regime of J < Q.25t. 
Therefore, they are neglected in what follows. 

In contrast, the spectral functions obtained with the 
usual spin t-J model show a dispersive feature at the 
highest energy, a similar weaker dispersive feature at low 
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FIG. 3: Ground state energy Eq as a function of the number 
of itinerant electrons for up to 5 holes (from top to 
bottom) added to a half-filled A'^ = 14 orbital chain, for Ejt = 
0, and for: (a) J = 0.25t, and (b) J = 0.125t. 



CL 



CL 




FIG. 4: Polarization of the orbital (filled circles, full lines) 
and spin (open circles, dashed lines) background P{n) Q at 
distance n from a single hole doped to a half-filled A'^ = 14 
chain, as obtained for the ground state of the orbital/spin t-J 
model without staggered fields (-Ejt = 0, hs =0), and for: 
(a) J = 0.25t {N^ = 7), and (b) J = 0.125t {N^ = 8). 



C. Correlation functions at finite doping 



energy, and incoherent spectral intensity between them 
[Fig. EJb)]. The two dispersive features have a period- 
icity of TT instead of 27r obtained for free electrons, are 
broadened and rather incoherent. Unlike in the 2D case 
where a hole is trapped and can move only by its cou- 
pling to quantum fluctuations in the spin background, 
leading to a quasiparticle behavior with a dispersion on 
the energy scale cx J,'"' the string potential is absent in a 
ID chain. The ID model shows spin-charge separation, 
because a hole may propagate after creating and leav- 
ing behind a single solitonic defect with two spins of the 
same direction next to each other 4ii^24^ However, the 
chain used in our case is too small for the spin-charge 
separation to be clearly visible. A more coherent compo- 
nent can be induced by introducing the coupling between 
the moving charge and spin fluctuations, which occurs in 
the presence of an external staggered field generating the 
string potentiali22i^ We address this issue in Sec. IIII El 

The present = 14 spin chain gives at half filling 
{S^S^_^_i) = —0.1491, a value being already quite close 
to the exact result —0.1477 for an infinite Heisenberg 
chain. Therefore, we expect that the results presented 
for the more classical orbital chain with partly localized 
wavefunctions are at least of the same quality, and suffer 
even less from finite size effects. 



First we consider the total energy Eq = (Hq) as a 
function of the number of \z) electrons at each doping 
X, and determine the actual distribution of Cg electrons 
in the ground state. The density of Cg electrons in two 
orbitals, rtj, and n^, depends on the ratio J/t. At half 
filling [x — 0) one finds — = 0.5, and the occupied 
orbitals alternate. At large J/t ^ 1 one expects that the 
orbital order is still close to perfect in the regions which 
separate orbital polarons, with the densities Uz — 0.5 and 
Ux — 0.5 — X. On the contrary, in the limit of J/t — > 
a single hole (in a finite chain) suffices to destabilize the 
orbital order, causing a transition to an itinerant (FO) 
state with Uz — 1 — x, and = 0. 

The energy obtained for an A^ = 14 chain filled by 
9 < A^ < 14 electrons for two values of J = 0.25t and 
0.125i are shown in Fig. |2| The tendency towards elec- 
tron delocalization is quite distinct already for a higher 
value of J = 0.25t, with the minimum of Eq moving to 
Nz > 7 with increasing doping x. The number of \z) 
electrons increases by one (to A^^ = 8) at the electron 
fiUing of A^e = 12 and 11 electrons. At A^e = 10 there 
remains just a single \x) electron which still blocks the 
hopping along the chain, while at A'e = 9 all electrons 
are in \z) orbitals, and one obtains a metallic state. This 
transition to a metallic state is faster at a lower value 
of J = 0.125i [see also Fig. ^c)]. Here adding a hole 
increases simultaneously the number of \z) electrons by 
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FIG. 5: Orbital correlations R{n) ^ at distance n from a 
single hole doped to a half-filled A'^ = 14 chain (filled circles, 
full lines), as obtained for the orbital t-J model for Ejt = 0, 
and: (a) J — 0.25t, and (b) J = 0.125t. Spin correlations ob- 
tained in the ID t-J model {hq = 0) is shown for comparison 
by open circles and dashed lines. 



one, and one finds — 8 and iV^ = 9 for iVg = 13 and 
Ne = 12, respectively, while doping by three holes gives 
already a metallic state with A^^ = A^e = H- 

The correlation functions in the ground state obtained 
after doping the chain by a single hole are very trans- 
parent at a higher value of J = 0.25t, and are easily 
accessible by looking at the correlation functions l|HJl~ 
(|10|l . As shown schematically in Fig. da), a hole replac- 
ing \x) electron is then confined to a three-site cluster, 
which can be deduced from the hole-orbital correlation 
function P{n) [see Fig. EJa)]. First of all, at the nearest- 
neighbor of a hole one finds preferably a \z) electron, with 
P{1) < —0.25 indicating that a hole spends more time 
at a central site than at either outer site of the three-site 
cluster [P(l) = —0.25 would correspond to the bonding 
state at J = and P(l) = —0.5 to a static hole at a cen- 
tral site]. This also causes weak alternation of Pin) with 
increasing n for further neighbors {n >2). This result is 
different again from the spin t-J model, where almost no 
preference for the spin direction is found already at the 
nearest neighbor of a doped hole. Here the oscillations 
between even and odd neighbors are also much weaker 
as a hole is now delocalized, and one is averaging over 
several different configurations when P(n) is evaluated. 
They correspond to the domains with opposite spins, and 
weak oscillations result here only from holon-spinon cor- 
relations. 

The ground state obtained for one hole at J = 0.125t 
is qualitatively different, as doping by one hole generates 



FIG. 6: Orbital correlations R{n) @ in a doped A*' = 14 site 
chain (filled circles, full lines), as obtained for the orbital t-J 
model at J = 0.125f, -Ejt = for doping by: (a) two holes 
(insulating phase), and (b) three holes (metallic phase). Spin 
correlations as in Fig. |3 



a single orbital flip and the number of itinerant \z) elec- 
trons increases to Nz ~ 8. Fig. El^b) shows that this 
defect in the otherwise perfect orbital order occurs close 
to the hole, increasing the island over which the hole 
can delocalize, extending now over five sites. Hence, one 
finds P(l) ~ -0.44 and P{2) ~ -0.25, and only starting 
from the fourth neighbor this correlation is weakly pos- 
itive. In contrast, decreased value of J does not cause 
any significant change in the P{n) correlations for the 
spin t-J model — they are only slightly weaker than for 
J = 0.25^. 

The orbital order, measured by second correlation 
function R{n) @, remains perfect at a sufficient distance 
from the doped hole, with R{n) = —0.25 [Fig. Efa)]. In 
case of a higher value of J = 0.25t, the order is perfect 
starting from n = 2, while for J — 0.125t with a larger is- 
land of the itinerant phase, it starts only at n = 4. In ad- 
dition, the orbital correlation i?(l) is here positive [Fig. 
EJb)], as two \z) orbitals occur frequently next to each 
other both at the first and at second neighbor of the hole 
which polarizes a larger five-site cluster. This behavior 
shows that phase separation occurs here in the regime 
of small J/t, and its mechanism which originates from 
orbital physics is completely different from that known 
from the spin t-J model at large J/t > 1. 

The spin t-J model has practically the same spin-spin 
correlations for both values of J. Due to quantum fluctu- 
ations in the Heisenberg ID chain, the correlation func- 
tion R{n) ~ —0.14 for n > 1 [Fig. EJa)] is much reduced 
from the classical value R{n) — —0.25 found for the or- 
bital model. In fact, the long-range order is absent in the 
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FIG. 7: Orbital and spin order (filled and empty circles) 
around a hole Z llOII as obtained for an A'' = 14 chain with 
increasing J/t doped by: (a) one hole, and (b) two holes. 
Other parameters as in Fig. |21 



ID spin chain, and the quantum fluctuations contribute 
to the energy (per site), which is much lower than in 
the orbital case, eo = J{S^ ■ = 3J(5f5,f+i). The 

present = 14 spin chain gives at doping by one hole 
(ni5'f^„5,f_^„^;^) = —0.1448 for n = 6, 7, a value which 
is only little reduced from —0.1491 found at half filling, 
showing that in spite of certain hole delocalization, spin 
intersite correlations are almost undisturbed by a single 
doped hole when the distance from it is sufficient. 

When doping increases, a gradual crossover in R{n) 
correlations towards a metallic chain with \z) orbitals 
occupied is found. The region doped by holes extends 
for two holes over nine sites and gives R{n) < start- 
ing from n — 3 [Fig. EJa)]. The correlation function 
R{n) is averaged over several bonds, and the values of 
R{n) ~ —0.17 for n = 6,7 result from a superposition 
of the classical order {{T^T^j^^) — —0.25) in the insulat- 
ing region and a positive value of {TfTl^^^) correlations 
within the metallic cluster. When the insulating phase 
disappears at doping by three holes [Fig. Elb)], R{n) > 
shows that only \z) orbitals are occupied. Three holes in 
the metallic phase avoid each other, and one finds a weak 
local minimum at n = 4 in R{n) correlations. The spin 
system gives instead negative R{n), as here the hole hop- 
ping does not destroy the AF spin order. The correlation 
functions R{n) for the spin t-J model have similar values 
for either two or three holes, while local minima found 
at n = 7 and n = A indicate again characteristic distance 
between the holes in both cases. 



The polaronic character of a doped hole is confirmed by 
the positive orbital correlation Z > between its neigh- 
bors 1)10(1 . At J < Q.Olt the ground state is metallic and 
Z = 0.25. Next, when J/t increases, the hole is first 
trapped in a large metallic cluster {Z ~ 0.22), with its 
size gradually decreasing down to a five-site cluster when 
J > O.OSt — then a hole is surrounded predominantly 
by two \z) electrons [Fig. [3a)]. This ground state, with 
Nz = 8 and = 5 electrons in the chain remains stable 
up to J = 0.20t. At J/t > 0.20 the ground state changes 
again to A^^ = 7 and Nx — 6 which gives a small orbital 
polaron [see Fig. ^a)]. Due to finite superexchange en- 
ergy J, the probability that the hole is in the center of 
a three-site cluster, with two neighboring \z) electrons, 
is somewhat higher than that it occupies either of side 
atoms, with one \z) and one \x) occupied orbital next 
to it. As a result, the correlation function Z ~ 0.02 is 
now weakly positive. For the filling by two holes one 
has again five different ground states [Fig. Efb)], with 
Z = 0.25, Z ~ 0.22, Z ~ 0.20, Z ~ 0.16, and Z ~ 0.07 
for increasing J/t. The transition to two separated three- 
site polarons occurs here at a somewhat higher value of 
J ^ 0.30t than the transition to a single three-site po- 
laron for the doping by one hole. 

In contrast, the spin correlations between two hole 
neighbors are negative in the spin t-J model, both for 
one and for two holes (Fig. 0)i showing the solitonic 
character of charge defect in the spin chaini^ At low J/ 1 
these correlations are more pronounced, while increasing 
J/t leads to more delocalized holes and to weaker corre- 
lations Z ~ —0.10 and Z ~ —0.12 for the states doped 
either by one or by two holes at J/t = 0.1. Somewhat 
stronger AF correlations in two-hole case confirm the soli- 
tonic character of charge defects in a ID spin chain, with 
a soliton compensated by an antisoliton. 



D. Spectral functions at finite doping 

Next we analyze the spectral functions obtained at fi- 
nite doping. Starting from the ground state of an orbital 
chain with a single hole, we found that in the region of 
very small J [J ^ Q.Olt for = 14 chain) the spectral 
function of the orbital model shows a free propagation of 
a hole within the entire chain, and electrons fill \z) or- 
bitals. This case (not shown) is however only of theoreti- 
cal interest, and at J ~ 0.125<, adequate for for LaMnOa, 
one finds instead fairly localized spectra [Fig. |HIa)]. One 
recognizes the maxima which correspond to localized \x) 
excitations at w — /i ~ — l.Si, and the structures corre- 
sponding to the bonding and antibonding states of the 
\x) excitations at w — /z ~ — 3.3t and — 0.3t. The spec- 
tra are fc-dependent and the spectral weight moves to 
higher energies with increasing momentum fc, following 
the tight-binding dispersion at E^'y = 0, £fe_ = — 2icosfc 
(|16|l . The total widths of the spectrum in Fig. |HIa) is 
close to At, i.e., to full free-electron dispersion obtained 
for electrons in \z) orbitals. 
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FIG. 8: (color online) Spectral functions A{k, ui) for = 14 
chains doped with one hole for J/t — 0.125: (a) orbital t- 
J model O with Nz = 8 and Nx =5 — solid and dashed 
lines for \z) and |a;) excitations; (b) spin t-J model 11211 with 
Ni — 6 and A'^x = 7 — solid and dashed lines for hole and 
electron excitations. Insets show the densities of states n(aj). 
Other parameters and peak broadening as in Fig. |5| 



The spectrum obtained for spin t-J model is drastically 
different. A broadened quasiparticle band crossing the 
Fermi energy resembles a tight-binding dispersion with 
a reduced bandwidth [Fig. |HIb)]. For low momenta k, 
where a larger system would show spin-charge separation, 
the spectra are especially incoherent. Most of the spec- 
tral intensity is in the incoherent part, but one recovers a 
coherent propagation in the limit of J ^ 0, which occurs 
in this case without polarizing the chain, as both spin 
flavors are mobile in this limit and the spin direction is 
irrelevant. 

As we have shown in Fig. Otb), doping = 14 chain by 
three holes is sufficient for the crossover to the metallic 
state for a realistic value of J/t — 1/8. The spectral 
function A{k,uj) confirms the metallic behavior of the 
orbital chain for this filling, and the excitations in |a;) 
orbitals contribute with a /c-independent maximum only 
above the Fermi energy [Fig. Ela)]. In contrast, the 
spectral functions in the spin t-J model [Fig. I^b)] are 
qualitatively very similar to those found at lower doping 



FIG. 9: (color online) Spectral functions A{k, to) for N = 14 
chains doped with three holes for J/t = 0.125: (a) orbital t-J 
model CJ with iV^ = 11 and Nx = 0, (b) spin t-J model 
with Ni — 5 and A^^ — 6. The meaning of solid (dashed) 
lines, insets and other parameters are the same as in Fig. |H| 



by one hole [Fig. |HIb)] , but with more weight transferred 
now to the quasiparticle states above the Fermi energy 
/i at A: > tt/2. The weight of the quasiparticle states for 
Lu > fi is enhanced, while the incohorent spectral weight 
is decreased for these values of k. 

The ground state of the orbital chain is insulating at 
a larger value of J = 0.25t [Fig. Ol^a)]. Therefore, the 
spectral functions consist in this case of several incoher- 
ent features for each k value, with the first moment of 
the spectra following again the single-particle dispersion 
(Fig. I10|l . The incoherent feature at w ~ /i, originating 
from excitations within \x) orbitals, is partly below the 
Fermi energy. This part of the spectral weight follows 
from |a;)-hole excitations which are still possible at the 
filling by A^^ = 9 and — 2 electrons. 



E. Polarons induced by the JT effect 

The energies of insulating and metallic phase are close 
to each other for typical parameters. Therefore, even a 
moderate JT energy Ejt = 0.25t is sufficient to stabilize 
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FIG. 10: (color online) Spectral functions A{k, cj) for an 
A'^ = 14 orbital chain doped with three holes for J/t = 0.25. 
Solid and dashed lines for \z) and \x) excitations. Inset, other 
parameters and peak broadening as in Fig. |21 



the insulating phase in a broad regime of doping, with 
occupied |x) orbitals fragmenting the chain into smaller 
units. When a single hole is doped, the JT field has qual- 
itatively a similar effect as a larger value of J — a hole 
is then well localized within a three-site cluster, and the 
hole-orbital correlation P(n) exhibits more pronounced 
alternation, indicating robust orbital order [Fig. Illf a)]. 
as in Fig. [S^a). The orbital intersite correlations R{n) 
are perfect already beyond the second nearest neighbor 
of the hole [Fig. Efb)]. 

When more holes are added, the orbital order gradu- 
ally softens, but a clear tendency towards orbital alter- 
nation is observed even for high doping with four holes 
[Fig. [TTT bl]. Note that a second hole may be added at 
a fourth or sixth site away from the first one. This, to- 
gether with weak hole hopping within the three-site clus- 
ters, weakens both hole-orbital correlations P{n) and the 
orbital order R{n) at n = 4 and (less) at n = 6. For the 
present case of = 14 sites these latter correlations are 
reduced more for n > 4 than for n = 2. With a third 
hole added, we found separated holes in three-site clus- 
ters, with each hole occupying predominantly the central 
site of its cluster in order to minimize the JT energy. A 
typical inter-hole distance is then four lattice constants, 
and for this reason P(4) becomes negative, and i?(4) is 
reduced stronger than R{n) for any other n > 0. The 
fourth hole is added at one of \x) orbitals within a longer 
\z) — \x) — \z) — \x) — \z) unit, and makes then the first 
bigger cluster, with two holes and three mobile \z) elec- 
trons. Particularly in this regime of parameters and up 
to this doping regime, the distribution of holes is remi- 
niscent of doping CuOa chains in YBa2Cu306+a;, where 
holes are doped first in separated Cu02 units, where they 
are trapped, generating fragmented units of CuOs chains 
and causing jumps and plateaus for the hole counts in the 
Cu02 planes as a function of dopingiMi 
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FIG. 11: Evolution of orbital correlations with increasing 
hole doping, as obtained for J — 0.125t and Ejt = 0.25t: 
(a) hole-orbital correlations P{n), and (b) hole-orbital-orbital 
correlations R{n). Filled circles, empty circles, filled triangles 
and empty triangles for one, two, three, and four holes doped 
to an A'' = 14 orbital chain. 



The spectral functions A{k,uj), presented for two rep- 
resentative doping levels with one and three holes in Fig. 
1121 show that they are remarkably similar in the entire 
regime of doping x < 0.3. Excitations in ja;) orbitals lead 
to maxima below and above the Fermi energy, with the 
spectral weight moving gradually to higher energies un- 
der increasing doping. For hole excitations in \z) orbitals 
one finds three nondispersive features: the central peak 
for doping within the orbital ordered regions, and two 
side peaks at energies corresponding to the excitations 
of orbital polarons. The spectral weight of the latter 
features increases on the cost of the central peak when 
doping increases, and the spectral weight distribution de- 
pends on momentum k. At doping increasing up to five 
holes (not shown) the spectra are still fairly localized, 
and the ja;) spectral weight is further reduced for lu < fi. 

A remarkable change of the spectral functions A{k,uj) 
has been found for the spin t-J model in presence of a fi- 
nite staggered field hg — 2 J. The field generates a string 
potential, and simulates thus a 2D model. ^""^ However, 
in this quantum model, in contrast to the orbital model 
shown in Fig. I12r a'l. a hole is not confined within a po- 
laron but can move coherently when it couples to quan- 
tum spin fiuctuations, which repair the defects in the spin 
background generated by a moving hole. As a result, 
a quasiparticle peak emerges close to the Fermi energy 
(Fig. [T^ , indicating a propagation of a hole dressed by 
spin excitations, with a dispersion ^ 4 J [note the units in 
Eq. (|12|) ]. The remaining spectral weight is distributed 
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FIG. 12: (color online) Spectral functions A{k, uj) for an insu- 
lating iV = 14 orbital chain for J/t = 0.125 and Ejt = 0.25t, 
obtained for increasing doping by: (a) one hole, and (b) three 
holes. Solid and dashed lines for \z) and \x) excitations. Insets 
and peak broadening as in Fig. |5| 



over incoherent features at lower energy while the low 
energy quasiparticle band visible in Fig. |S1 has vanished 
as well as the incoherence of the band for small k. 

Moreover, the quasiparticle band has lost a simple 
tight-bindinglike dispersion found before at hg = 0, and 
instead of crossing the Fermi energy, it has periodicity 
TT because of the doubled unit cell and folds back to- 
ward lower energies after reaching the Fermi energy at 
k = 7r/2. The peaks are sharpest near this point, con- 
firming the quasiparticle character of this excitation, but 
have lower spectral weight at fc = tt than at fc = 0. The 
quasipartiple band close to the Fermi energy leads to a 
quite distinct peak in the density of states (see inset of 
Fig. 1131) . It is separated by a pseudogap from the elec- 
tronic excitations at w > /i, and also by another pseudo- 
gap at ~ 4 J below the Fermi energy from the incoherent 
part of the spectrum at lower energies. Thus, the present 
case is radically different from the case of /is = [see 
Fig. |2Ib)], and resembles the spectral functions for a 2D 
modelM 
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FIG. 13: (color online) Spectral functions A{k,uj) for an = 
14 spin t-J chain p2|l at doping by one hole, as obtained for 
J/t = 0.125 with a staggered field hs — 2J . Solid and dashed 
lines indicate hole and electron excitations. Inset shows the 
density of states n(tj). 



IV. ORBITAL POLARONS AT FINITE 
TEMPERATURE 

A. Effective orbital t-J model 

At increasing temperature the FM spin order assumed 
in Sec. EJis gradually destroyed and AF spin configura- 
tions on the bonds occur with finite probability, modify- 
ing the form of both i and J terms in Eq. 1^. Even when 
the ground state at T = is FM, superexchange interac- 
tions which originate from t2g electron excitations play 
an important role and contribute at finite temperature. 
These interactions are frequently treated as an effective 
AF superexchange between core S* = 3/2 spins, but de 
facto they depend on the total number of d electrons of 
two interacting Mn ions^ We have verified, however, 
that the t2g superexchange terms derived for these dif- 
ferent configurations are of the same order of magnitude, 
so it is acceptable to consider their effect as equivalent to 
a Heisenberg interaction with an average exchange con- 
stant J' > 0. Therefore, we include in the present ID 
model the spin interaction. 



H.J' — J' ^ [Si ■ Si 



+1 



(18) 



The superexchange due to Cg electron excitations con- 
tains spin scalar products multiplied by orbital inter- 
actions on the bonds, and the full many-body problem 
would require treating the coupled spin and orbital dy- 
namics. Here we decouple spins and orbitals in the mean- 
field approximation, and study the orbital correlations 
and their consequences for the magnetic order by replac- 
ing the scalar products of spin operators on each bond 
by their average valuesji*^ 



1 



(19) 
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where = cos(0i^i+i/2)e*''^'''+'^ , and directions of two 

classical spins at sites i and i + 1 differ by angle ^i^^+i. 
The complex phase Xi,i+i does not have any effect in the 
present ID model<i£ 

We investigate the effective orbital t-J model, 



niS) = Ht + Hj + Hj,+HjT, 



(20) 



where the AF interactions of core spins H18() and the JT 
field acting on orbital variables (O are playing the role 
of external fields. The hopping term Ht describes the 
dynamics of Cg electrons which are locally coupled to t2g 
core spins by large Hund's exchange Jh element, and 
have their spins parallel to them in the ground state. 
In agreement with double exchange mechanism,^'™ the 
hopping along each bond {i,i + l) is then modified by the 
spin order to 



tu. 



(21) 



and vanishes when the spins are antiparallel. This ap- 
proximation allowed to establish the existence of mag- 
netic polarons in the ID Kondo modeliii 

The superexchange due to Cg electrons cx J gener- 
ates orbital interactions in Hj which follow from vir- 
tual charge excitations, either for Mn'^+-Mn'^+, or for 
Mn^+-Mn'*+ pairs. Consider first the spin-orbital model 
of Ref. ^2d for an undoped LaMnOa. So far, we included 
the orbital superexchange [Eq. Q for a situation when 
high-spin eg excitations dominate, but low-spin excita- 
tions would also contribute at finite temperature. Tak- 
ing realistic parameters for the Coulomb interaction U 
and Hund's exchange Jh, the low-spin '^Ai, '^E, and "142 
states have the energies: 7.3, 7.8, and 9.6 eViSi We keep 
the relative importance of various excitations and write 
the terms resulting from various charge excitations to the 
superexchange, after averaging them over the spin con- 
figuration, as follows: 
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(22) 



In a doped system one finds only AF superexchange 
terms for Mn3+-Mn'*+ bonds [while FM ones are explic- 
itly included as double exchange in the orbital model 
()20|l ] , which contribute only if an Cg electron at site Mn'^+ 
occupies a directional orbital along the considered bond 
direction, in our case a 3z^ — orbital, and can hop to 
its Mn^+ neighbor. Taking the excitation energy to the 
low-spin configuration e{'^E) ~ 3 eVir^ one finds 
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{Si ■ Si+i — 3) 



1). (23) 



Using the coefficients given in Eqs. (|22|l and (|23|l . one 
arrives at the effective orbital t-J model (jJOJ, with 



= J^E(2"L+i + 3)(2i;^i;^+i-in.n.+i 



Uis+ijnizTii+i^^^ 



(24) 



(25) 
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The hopping {Ht) and the interaction {Hj) terms depend 
on the actual spin configuration which fixes the bond 
variables {ui^i+i} [see Eq. (I19|) ]. In the FM state at 
T = all Ui^i+i — 1, the AF terms vanish, and one 
recovers the form of Ht and Hj used in Sec. 

It is a crucial feature of the effective orbital model 
given by Eq. H20() that spin interactions are influenced 
by orbital correlations along the chain, and the latter can 
support either FM (for alternating x/z orbitals) or AF 
(for polarized z orbitals) spin order. At T = one can 
find on optimal state by minimizing the total (internal) 
energy of the system, E = {H{S)), over the spin and or- 
bital configurations. In fact, at T = two solutions are 
possible, depending on the parameters in Eq. H20|l . Let 
us consider first a purely electronic model with Ejt = 
at half filling (we will show below that the situation is 
similar at Ejt > 0, but the region of stability of the FM 
states is extended). If J' = 0, FM order is stable and 
coexists with alternating orbital order. This situation 
was discussed in detail in Sec. IIIII However, a relatively 
small value of J' = 0.0125t is sufficient to compensate 
the energy difference between this state and an AF state 
with occupied \z) orbitals, and (at T = 0) one finds the 
latter state for J' > 0.0125t. It is interesting to inves- 
tigate a competition between these magnetic states at 
finite temperature when the chain is doped. 

We investigated the spin and orbital correlations for 
the effective orbital t-J model H20|l at finite temperatures 
by employing a combination of a Markov chain Monte 
Carlo (MCMC) algorithm for the core spins, with Lanc- 
zos diagonalization for the many-body problem posed by 
the orbital chain of = 12 sites with periodic boundary 
conditions, at each given distribution of classical vari- 
ables {iti_i+i} for i — 1, • • • , 12, which stands for a par- 
ticular spin configuration. The partition function to be 
then evaluated is 



Z = / V[S] Tic e 



-/3W(S) 



(26) 



where /? — l/ksT (we adopt the units with fc^ = 1), 
and / 2? [5] denotes the integral over the A^-dimensional 
space of all core spin configurations for the chain of 
length N. The core spins S = {Si} determine the 
site-dependent hopping parameters {ui^i+i} and thus the 
fermionic Hamiltonian 7i(iS) is fixed. 



The trace over the fermionic degrees of freedom, 
Tr^e"''^^'^^ =: w{S), gives the statistical weight for S 
and is sampled by the MCMC. For independent electrons 
it can easily be evaluated;^'^ for interacting electrons one 
has to use a Lanczos algorithm/^ but as w{S) is strictly 
positive, one still has no sign problem. Since the posi- 
tion of the \x) electrons is conserved, the Hamiltonian 
(|20|l has a block-diagonal structure, with each block cor- 
responding to one fixed distribution of |a;) electrons along 
the chain. There are, of course, many such blocks, but 
their subsequent diagonalization is still much faster than 
the diagonalization of the complete matrix. For the ex- 
treme case of the completely filled chain, each block has 
dimension one, meaning that the Hamiltonian is already 
diagonal and no matrix-vector multiplication has to be 
performed. Furthermore, degenerate eigenvalues can be 
resolved, if they are in different blocks. 

For the MCMC updates, 'w{S) was calculated from the 
lowest few eigenenergies of each block for an TV = 12 
chain, until a new Lanczos step did no longer modify the 
contribution from this block. Observables were only cal- 
culated from the lowest 14 eigenstates of the whole space. 
In order to monitor this approximation, the Boltzmann 
factor of these states was measured. The weight of the 
highest included state was approximately 1.5 percent for 
the worst case (filled chain, pt = 20, J' = 0.02, Ejt = 0), 
0.3 — 0.5 percent for f3t = 50 and the filled chain and neg- 
ligible for /3 = 100 or finite doping. This means that for 
those observables, which are calculated from the eigen- 
states of the Hamiltonian, at most a few percent of the 
total weight were missed. For the MCMC updates and 
for observables which do not need the eigenvectors (core 
spin correlation and total number of |a:)/|z) electrons), 
the error was even smaller. 

The core spins were rotated for whole sections of the 
chain at once. Because acceptance for core spin updates 
was high, we performed two or three such rotations be- 
fore testing acceptance. Every spin was therefore rotated 
several times per sweep. The numbers and Nz of 
\x) and \z) electrons remained fixed for the evaluation 
of 'w{S), and every five updates we proposed to increase 
(decrease) and decrease (increase) N^, respectively, 
thereby sampling and . The total number of elec- 
trons Ne was kept fixed. Between measurements, 40 to 
100 lattice sweeps (depending on temperature) were done 
in order to decorrelate the samples. We then employed 
autocorrelation analysis and found the samples to be un- 
correlated. To reduce the statistical errors, 200 uncor- 
related samples were obtained for each set of parameter 
values. 



B. Orbital order versus spin order 

In addition to the correlation functions studied in Sec. 
mil spin correlations (S^Sj) are now investigated by eval- 
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FIG. 14: Orbital correlations (TfTf+i) ^ at half filling, as 
obtained at low {pt — 100, open circles), and at intermediate 
(/3t = 50, filled circles) temperature for an A'' = 12 site orbital 
chain with J = 0.125t, and: (a) J' = Ejt = 0, (b) J' = 
0.02t, Ejt = 0, and (c) J' = 0.02t, Ejt = 0.25t. Statistical 
errors are smaller than the symbol sizes. 

uating the spin structure factor, 

which depends on the ID momentum fc, and follows from 
Monte-Carlo simulations of = 12 chains. The averages 
{{S^S^^i)} are calculated from the spin configurations 

S = {Si} determined by the MCMC updates for {ui^i±i}. 

The interplay between spin and orbital correlations be- 
comes transparent by varying the core spin AF superex- 
change J' ()18|l due to t2g electrons. We investigated spin 
and orbital order for two characteristic values of temper- 
ature: (3t = 100 and (3t = 50, corresponding to T 60 
and ~ 120 K ioi t ^ 0.5 eV, i.e., well below the magnetic 
transition. First we consider the case of J' = which 
reproduces the ground state analyzed in Sees. ITlland lTTll 
While the orbital alternation, measured by T{n) corre- 
lation function lO, is perfect at T = 0, it softens some- 
what when temperature increases, but is still robust at 
temperature f3t — 50, as shown in Fig. EJa). Orbital 
alternation supports the FM spin order at half filling, 
which gives a distinct maximum at fc = of the spin 
structure factor S{k) (see Fig. I15II . In the weak doping 
regime with one or two holes added, the spin correla- 
tions, (niiS'?^„S'f_^„_,_j), are driven by superexchange and 
are FM at any distance n from the hole, and only weakly 
depend on n. These correlations increase by a factor 
close to two when the doping changes from two to three 
holes. This explains why the maximum of S{k) at fc = 
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FIG. 15; (color online) Spin structure factor S{k) as obtained 
for an undoped A'^ = 12 orbital chain Il2fl|l . and for increasing 
hole doping up to five holes: (a) at low temperature /3 = 
100, and (b) at intermediate temperature /3 = 50 (in units 
of Statistical errors are smaller than the symbol sizes. 

Parameters: J = 0.125t, J' = 0, Ejt = 0. 



remains almost unchanged in the low doping regime by 
one or two holes, but is next strongly enhanced when 
doping increases to three holes. Precisely at this concen- 
tration electrons redistribute within the chain and oc- 
cupy practically only \z) orbitals, giving a metallic state. 
This demonstrates that double exchange plays a primary 
role in the observed insulator-metal transition and sig- 
nificantly enhances the stability of the FM order in the 
metallic phase. 

At J' = 0.02t the orbital correlations found in the half- 
filled chain are markedly different [Fig. I14f b)]. At low 
temperature {f3t = 100) they indicate that primarily (but 
not only) \z) orbitals are occupied, while no \x) electrons 
were found in the ground state at T — 0. This orbital 
state is induced by finite J', and supports the AF spin 
interactions due to Cg excitations which select then low- 
spin states in the Eg excitations along \z) — \z) bonds. 
When a single hole is doped, the spin correlations remain 
still AF at low temperature [Fig. llBf a)]. except for the 
spin-spin correlation between the site occupied by the 
hole and its nearest neighbor, giving rise to a small FM 
polaron. Therefore, the value of S{0) increases already at 
this low doping, and the electron density in \z) orbitals 
is enhanced close to the hole itself. 

When temperature increases, it becomes clear that the 
orbital |a;)/|2;) alternation, supporting the FM spin inter- 
actions, competes at half filling with the above uniformly 
polarized chain with occupied \z) orbitals, supporting the 



FIG. 16: (color online) Spin structure factor S{k) for an TV : 
12 orbital chain ^ as in Fig. ^ but for J' = 0.02t. 



AF spin order. A clear tendency towards orbital alter- 
nation is detected by a negative nearest-neighbor orbital 
correlation T(l) ~ -0.10 [Fig. Elb)] at higher tem- 
perature {f3t = 50). Although the value of S'(7r) is still 
larger than S{0) at half fiUing, both (weak) maxima, cor- 
responding to FM and AF order, become almost equal 
already for doping with one hole [Fig. I16f b)] . This case 
is qualitatively similar to a chain doped by two holes at 
(3t = 100 [Fig. llGf a)]. where also two (stronger) max- 
ima of S{k) indicate coexisting islands of FM and AF 
spin correlations along the chain. Thus, we found that 
the AF correlations are gradually changed into FM ones 
with increasing temperature. This trend follows from the 
difference in the energy scales — when the thermal mag- 
netic excitations destroy the energy gains due to J and 
J', the kinetic energy oc t is a dominating energy which 
can be optimized by selecting \z) occupied orbitals and 
FM spin correlations. 

At doping by three holes the spin correlations are 
FM, both at low {/3t — 100) and at higher temperature 
{Pt ~ 50), showing that the double exchange enhances 
the effective FM interactions and changes the character- 
istic temperature at which the spin correlations weaken 
— thus the Curie temperature Tc would increase in a 3D 
case. This agrees with the experimental observations — 
indeed, the Curie temperature increases with hole doping 
X in the metallic regime.^^-^^ 

Although at T = one finds indeed the AF ground 
state with = 12 for a half-filled TV = 12 chain, 
many excited states with a few [a;) orbitals occupied are 
found at low energy and contribute already at f3t = 100. 
As a result, electrons are redistributed over Cg orbitals 
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TABLE L Total number of electrons in itinerant orbitals 
(top part) and in localized orbitals A'^^, (bottom part) for 
different doping level x and for increasing temperature fit, as 
obtained in the Monte-Carlo simulations for an A^ = 12 orbital 
chain [Eq. I|20|l ] (statistical errors are also given). Parameters: 
J = 0.125t, J' = 0.02t, £jT = 0. 















a; = 


X = 1/12 


X = 1/6 


X = 1/4 


100 


9.97 ±0.09 


9.80 ± 0.06 


9.64 ±0.04 


9.0 ±0.0 


50 


7.85 ± 0.07 


8.51 ±0.06 


8.81 ±0.04 


8.98 ±0.01 


30 


7.47 ± 0.07 


8.05 ± 0.06 


8.53 ±0.05 


8.88 ±0.03 


20 


7.34 ± 0.07 


7.85 ± 0.06 


8.18 ±0.06 


8.55 ±0.04 






A, 






Pt 


a; = 


X = 1/12 


X = 1/6 


X = 1/4 


100 


2.03 ±0.09 


1.20 ±0.06 


0.36 ±0.04 


0.0 ±0.0 


50 


4.15 ±0.07 


2.49 ±0.06 


1.19 ±0.04 


0.02 ±0.01 


30 


4.53 ± 0.07 


2.95 ±0.06 


1.47 ±0.05 


0.12 ±0.03 


20 


4.66 ± 0.07 


3.15 ±0.06 


1.82 ±0.06 


0.45 ± 0.04 



by thermal excitations, and the orbital polarization at 
(it = 100 is far from complete, with Nz — 10 and 
Nx ~ 2 (Table P). When temperature increases further 
to (3t = 50 {(3t — 20), orbital disorder increases and one 
finds Nz ~ 7.8 and ~ 4.2 {N^ ~ 7.3 and ~ 4.7). 
This demonstrates that the balance between FM and AF 
terms realized at half filling is rather subtle — increas- 
ing temperature favors more disorder in the chain which 
destroys a uniform AF phase, supported by charge exci- 
tations along \z) — \z) bonds. 

The population of localized \x) states in the chain de- 
creases quite fast with doping, particularly in the regime 
of pt < 100. At a; = 1/4 one finds almost no \x) defects 
in the metallic chains, except at rather high temperature 
Pt — 20. Increasing electron density in \x) orbitals with 
increasing temperature may be seen as a precursor effect 
for the metal-insulator transition at Tc in the interme- 
diate doping regime. 



C. Polaronic features at finite doping 

The spectral functions A{k,uj) obtained for the FM 
phase at half filling for J' = (not shown) are similar to 
those discussed in Sec. IIII Dl Apart from some broad- 
ening due to finite T, the one-hole excitations are again 
fairly localized, similar to those shown in Fig. [5] In 
the doped regime one finds first the localized spectra at 
low doping, similar to those of Fig. ISfa), with a distint 
pseudogap at the Fermi energy fi. As found before at 
T = 0, doping by three holes suffices for the crossover to 
the metallic phase, with somewhat broadened peaks in 
A{k,uj), following the one-particle dispersion due to the 
hopping within \z) orbitals, like in Fig. IHIa). 

In contrast, the spectral properties obtained for the 
AF phase found at half filling with J' = 0.02t are quite 
different (Fig. I17|l . First of all, there are predominantly 
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FIG. 17: (color online) Spectral functions A{k,uj) for a half- 
filled A = 12 orbital chain at temperature /3t = 100. 
Solid and dashed lines for \z) and |a;) excitations. Inset shows 
the density of states n{u)). Parameters: J = 0.125t, J' = 
0.02t, EjT = 0. 



\z) electrons at low temperature f3t = 50, and still more 
so at (3t — 100, so the spectral weight of the localized 
|a;) hole excitations, found at energies uj ^ — 1.3i and ^ t 
[Fig. IHT b)]. appears to be very low. The \z) excita- 
tions have insulating character and are incoherent. They 
are given by a superposition of two features, namely a 
dispersionless peak at the hole binding energy uj —2J 
and a weakly dispersive band with bandwidth ~ 0.5t. 
The first feature stems from the few \z) electrons which 
are surrounded by \x) electrons and is similar to that 
discussed in Sec. IIII Bl The second feature comes from 
electrons moving within the \z) polarized parts, with the 
strongly renormalized hopping cx ti i^i H21|) for aligned 
\z) orbitals along an almost perfect AF bond. If the AF 
order were perfect, the bandwidth would vanish, because 
no hopping would then be possible. Both the dispersive 
and the dispersionless structures are broadened owing to 
thermal fluctuations of the core spins. The superposition 
of these two bands yields a single peak with large inten- 
sity for k = 0, because the bandwidth is approximately 
~ 0.5t = 4 J and the both bands therefore coincide at this 
point. For the same reason, the spectral weight at fc = tt 
consists of two almost symmetric maxima. This inter- 
pretation is corroborated by the results for pt = 50 (not 
shown), where (i) the dispersionless feature has a higher 
weight, because the population of orbitals increases 
with increasing temperature (see Table I), and (ii) the 
structures in A{k,Lu) are still more broadened because of 
the larger thermal fluctuations. 

The spectrum changes rapidly when doping increases. 
At doping of X = 1/12 (one hole in an A^ = 12 chain) a 
polaronic peak is found above ^, while below /i the spec- 
trum separates into a broad incoherent part at interme- 
diate energies, and a low-energy peak with large intensity 
for low values of fc = and fc = 7r/6, being a symmet- 
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FIG. 18: (color online) Spectral functions A{k, ui) for an = 
12 orbital chain II2UI I at temperature pt = 100 doped by: (a) 
two holes, and (b) three holes. Solid and dashed lines, inset 
and parameters as in Fig. 1171 



ric image of the polaronic peak (not shown) . As de- 
creases (Table I), the \x) intensity below /i drops in doped 
chains. At doping x — 1/6, |a:) electrons were found only 
in some samples at low temperature {/3t = 100). It seems 
that these fluctuations alone would not suffice to make 
the system insulating, but AF correlations are here still 
strong (Fig. I16|l and generate a well-formed polaronic 
peak in the density of states for the electronic exciations 
(with oj > separated by a pseudogap from the inco- 
herent spectrum due to \z) hole excitations [Fig. llHr a)]. 
In this case the (localized) |a:) excitations are predomi- 
nantly electronic, leading to a distinct nondispersive peak 
above /i. Except for these latter excitations, the situa- 
tion is here very similar to that found in the one-orbital 
modelfii where FM polarons, generated dy doping, were 
embedded into an AF background. 

For the present parameters, doping x — l/A (by three 
holes) gives an almost perfect metallic chain (occupied 
\z) orbitals) at low temperature f3t = 100, and the \z) 
excitations are then coherent [Fig. I18f b')]. One finds 
seven excitations obtained for different momenta which 
sum up to five distinct structures in the density of states 
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FIG. 19: (color online) Spin structure factor S{k) as in Fig. 
1151 but for up to four holes. Parameters: J = 0.125t, J' = 
0.02t, EjT = 0.25t. 



n{Lu) — the two side features consist of two peaks each, 
representing joined intensities for /c = 0, 7r/6 and for k = 
57r/6,7r, respectively. An apparent pseudogap at ui — fi 
is a finite size effect and would disappear at large system 
sizes. 



D. Jahn- Teller effect at finite temperature 

We have shown above that the orbital and magnetic 
order are interrelated and influence each other. There- 
fore, as pointed out before,^*' not only the superexchange, 
but also purely orbital interactions which follow from the 
JT effect are of importance for the observed magnetic A- 
AF order in LaMnOs. Now we will show that also in 
the present ID model the JT effect may modify the mag- 
netic order. This is particularly transparent by consider- 
ing a uniformly polarized insulating state with occupied 
\z) orbitals, coexisting with an AF order, and stabilized 
by weak superexchange J' = 0.02t (Fig. I16|l . If an alter- 
nating JT potential, given by Eq. (0, increases, the AF 
order is easily destabilized — already at Ejt — 0.25t we 
found an almost perfect orbital staggering in the broad 
temperature regime at half filling [Fig. I14r c)]. which in- 
duces instead the FM spin order, visible as a broad max- 
imum in S{k) centered at /c = (Fig. 119(1 . However, this 
maximum is less pronounced and other correlations S{k) 
with k > are also present, unlike for J' = Ejt = 0, 
showing that the FM order induced by the JT potential 
is definitely much weaker. 

Even for moderate JT potential Ejt — 0.25t, increas- 
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FIG. 20: (color online) Spectral functions A{k, ui) for an = 
12 orbital chain l|2()|l . as obtained at temperature f3t — 100 
for: (a) three holes, and (b) five holes. Solid/dashed lines 
and inset as in Fig. El Parameters: J = 0.125f, J' = 0.02t, 
i;jT = 0.25t. 



ing doping does not change the spectral properties quali- 
tatively, and they remain dominated by localized excita- 
tions. For doping with one hole (not shown), two local- 
ized I a;) excitations are accompanied by three localized 
\z) excitations: the central peak corresponding to a lo- 
calized hole between two \x) electrons, and two satellite 
structures, a hole excitation at energy w — /j, ~ —2t (with 
large intensity for low k values), and an electron excita- 
tion at energy lu — fj. 0.5t, mainly contributing for large 
k values. These structures are similar to those found be- 
fore at T = [Fig. ITW a)]. but they are now broadened 
by thermal spin fluctuations. When doping increases, the 
intensity of \x) hole excitations decreases and moves to 
the Fermi energy lu — ii, while a satellite peak grows for 
electron excitations (uj > fi), as shown for doping x — 1/4: 
in Fig. I2ur a). At this doping the FM correlations are en- 
hanced by a factor ~ 1.5 with respect to the undoped 
case, but are still much weaker than those found before 
at the same doping for Ejt = (Fig. I16|l . where the 
FM phase was metallic. This demonstrates that the in- 
tersite magnetic correlations and the energy responsible 



for them are weaker in insulating FM manganites than 
in the metallic ones at the same doping. A similar trend 
was indeed observed for the values of the Curie tempera- 
ture TcA and one expects that also spin stiffness should 
be reduced by the JT coupling in CMR manganites 4i 

The transition to metallic phase is damped by the JT 
distortions, and although the FM correlations increase 
significantly at doping x = 1/3 (Fig. I19p . the chain re- 
mains insulating. Only for as high doping as x — 5/12, 
the \x) electrons are practically eliminated in the con- 
sidered temperature range f3t > 50, and a metallic be- 
havior takes over. This metallic state gives almost free 
dispersion of \z) states in A{k,uj), while \x) (electron) 
excitations are again localized [Fig. EUT b)] . 



V. SUMMARY AND CONCLUSIONS 

The present study clarifies that orbital degrees of free- 
dom are of crucial importance for the understanding of 
magnetic correlations in CMR manganites. First of all, 
for the realistic parameters of manganites the FM and AF 
state are nearly degenerate at half filling. In both cases 
the decisive term stabilizing the magnetic order origi- 
nates from the eg superexchange. Although the ground 
state at T = would be FM in the absence of the AF 
superexchange between core spins, is is easy to flip the 
balance of magnetic (and orbital) interactions and stabi- 
lize instead the AF order in a purely electronic model. 
Here we adopted the AF interaction J' = 0.02t between 
core spins in order to stabilize the AF spin order in an un- 
doped chain (the true value of J' ~ O.OOAt in LaMnOa, es- 
timated from the value of Neel temperature of CaMnOa, 
is smaller by a factor close to five^P)^ and to demonstrate 
a gradual crossover from an AF insulator to a FM metal 
under increasing doping. In this regime of parameters 
the present ID chain stands for the AF order along c 
axis in the A-AF phase realized in LaMnOs. 

We have shown that even in the ID model the mag- 
netic interactions are internally frustrated, with compet- 
ing FM and AF terms in the superexchange. A delicate 
balance between these terms is easily disturbed by the 
JT potential originating from the lattice. This allows 
one to investigate within the same framework both types 
of magnetic order which coexist in the A-AF phase of un- 
doped manganites. On one hand, if the oxygen distortion 
have almost no influence, or if a uniform polarization in 
\z) orbitals would be induced by them, the AF spin order, 
found along c axis in ^-AF phase of LaMnOa, would be 
supported. On the other hand, if the oxygen distortions 
are alternating, they induce alternating orbital order, as 
it happens in (a, b) planes of LaMnOs, and the FM terms 
are selected from the Cg superexchange. Note, however, 
that the FM correlations are weakly reduced from their 
values found in the absence of the JT interations (J' = 
and EjT = 0), demonstrating that the FM interactions 
do not depend explicitly on the JT terms, but are only 
induced by a given type of orbital order. In this way the 
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JT coupling to the lattice helps to remove the frustration 
of magnetic interactions in CMR manganites. 

The evolution of spin correlations for increasing dop- 
ing showed that indeed two mechanism are responsible 
for ferromagnetism: when the weakly doped ID chain 
is insulating, the FM interactions are induced both by 
the superexchange terms following from the high-spin ex- 
cited states, and by the local double exchange within 
polaronic states around single holes trapped in the insu- 
lating phase. The double exchange interactions is much 
stronger than the superexchange, and it fully takes over 
and operates in the metallic phase at higher doping. The 
difference between these two mechanisms is reflected by 
a fast increase of FM correlations at the insulator-metal 
transition which was investigated both within a purely 
electronic model, and including the JT potential induced 
by the lattice. 

It is quite remarkable that orbital polarons found in 
the present model with orbital degeneracy in the regime 
of AF spin correlations resemble FM polarons which oc- 
cur in the Kondo model^ This provides some support to 
a simplified picture of a nondegenerate conduction band 
which is able to capture the essential physics when the 
orbitals are polarized, and the orbital degrees of free- 
dom are quenched and do not contribute in any signifi- 
cant way. A conservative point of view, based on double 
exchange mechanism, is that the FM polarons compete 
with the AF order and cause a transition to the metallic 



FM phase. Yet, this is not the only possibility — we have 
shown that the FM phase at low doping could be insu- 
lating due to immobile orbital polarons, which allow to 
understand why this phase could be FM and insulating 
at the same time. Such polarons are expected to play an 
essential role in the insulator-metal transition within the 
FM phase in manganites. 

We believe that many qualitative features found in the 
present ID study are generic for the interplay between 
orbital and magnetic order in CMR manganites. Work is 
in progress on higher dimensional systems. Among oth- 
ers, an interesting question is to what extent the orbital 
order is modified when two Cg orbitals start to fluctuate 
more strongly in either 2D or 3D systems, both due to 
quantum effects and due to increasing doping. 
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